Brownian motion: a case of temperature fluctuations 
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A diffusion process of a Brownian particle in a medium of temperature T is re-considered. We 

^— V ' assume that temperature of the medium fluctuates around its mean value. The velocity probability 

f^ , distribution is obtained. It is shown that the stationary state is not a thermodynamic equilib- 

fv^ ' rium state described by the Maxwell distribution. Instead a nonequilibrium state is produced by 

_, ' temperature fluctuations. 



^ 



0^ 
(N 






o 



> 

00 

o 



X 



I. INTRODUCTION 



Generalized statistics or 'superstatistics' occur in non-equilibrium systems as a result of parameter (temperature, 

friction, energy dissipation, pressure, chemical potential, etc.) fluctuations [i|. An example of superstatistics is 

the Tsallis statistics in nonextensive statistical mechanics .2j. One of a dynamical realization of this statistics has 

^ . been constructed by a Langevin equation for the Brownian particle |3, l^j with the inverse temperature being a 

fluctuating parameter. In the paper we consider a more natural model with fluctuating temperature instead of 

its inverse. Fluctuations of temperature can play a significant role in many processes and phenomena. E.g., in 

"t^ astrophysics, the spectrum of temperature fluctuations of the cosmic microwave background radiation can change our 

-j-J view on the universe at epochs from redshifts of the order of ten thousand to nearly the present and can provide 

. ' important clues to inflationary models and the dark matter-energy problem u\. In plasma physics, an experimental 

jrt \ evidence of substantial temperature fluctuations has been found in mechanisms responsible for anomalous transport 

Cj . observed in tokamaks and stellarators Q- The concept of temperature fluctuations is used in the theory of heavy 

I ■ ion collisions and multiparticle production \(\. In the Rayleigh-Benard convection, temperature fluctuations can be 

'^ ■ passively transported in the turbulence regimes Q- Characteristics of temperature fluctuations in living tissue has 

^ ■ been studied in Q . Below, we study the influence of temperature fluctuations on motion of the Brownian particle. 

O , As mentioned, the similar problem has been studied previously, mainly in the context of the Tsallis statistics |2j 

with application to velocity fluctuations in a turbulent flow [J. However, the studies have been limited to inverse 

temperature fluctuations and to the 'static' case when fluctuations are represented by a time- independent random 

variable. The problem is mathematically trivial in the sense that the probability density of velocity can be calculated 

for an arbitrary random variable modeling fluctuations. The more realistic model seems to be the 'dynamic' model 

based on time-dependent noise for which temperature fluctuations are represented by a stationary stochastic process. 

Then, as is shown below, the problem becomes non-trivial, even for the simplest model of temperature fluctuations. 

\^ , Let us remind that in the classical theory of diffusion, a position x = x{t) of a one-dimensional motion of a 

f^ ' Brownian particle of mass m moving in an equilibrium homogeneous medium of temperature T is described by a 

TJ" ] Newton equation with a random force which, according to the fluctuation-dissipation theorem, has the form [l^ 

O ' 



mx + 'yx^y/2-fkT£,{t), (I) 



where 7 is the friction coefficient (given by e.g. the Stokes formula), k is the Boltzmann constant and ^(t) is a random 
I ' force modeled by the Gaussian white noise, 
^ ■ 

g: <^(i)>=o, <mas)>^s{t-s). (2) 

^ ' The velocity u = i is the Ornstein-Uhlenbeck stochastic process and its probability density P{v,t) obeys the Fokker- 
^ Planck equation 



dPiv,t) _ 7 dP{v,t) jkTd^P{v,t) 
dt m dv m? dv"^ 



^2 '^'^ "^"^ ^^ ' ^"^^ Wi 1 ^ 



A general solution of this equation is given by the expression 

/oc 
p(i;,i|«o,0)F(wo,0)dz;o, (4) 

-00 

where P(i;,0) is an initial distribution and the transition probabihty distribution 

p{v,t\vo,0) ^ [2na\t)Y^'\^^\~\v - voe'-^"^Y I2a^{t)\ . (5) 



The variance 

kT 



CT2(t) = — fl-e-27*Any (6) 

m \ J 

The stationary velocity distribution function Pm{v) does not depend on the initial distribution P(v,0) and is the 
Maxwell distribution, 



It means that the stationary state in the velocity space is a thermodynamic equilibrium state. 

Now, let us consider the situation when the fluid is in a nonequilibrium steady state. What is a stationary state of 
the Brownian particle in this case? We should use a theoretical framework which enables to answer this question in 
a unified manner. For the moment, however, such a universal theory does not exist. First of all, we should specify a 
nonequilibrium state of the fluid: We assume that this state is not far from equilibrium and can be described similarly 
as an equilibrium state with the only exception that now temperature T is time-dependent, T = T(t). We follow and 
extend the proposal of Beck [ll Q from adiabatic to non-adiabatic temperature changes and app ly Eq. (^ to this 
case. The problem is whether and when we can use Eq. Q if T = Tjf), cf. also polemics in |lll |. Let us remember 
that when (Q is derived from the microscopic Hamiltonian model jl2| . it is assumed that the fluid (medium) is in 
the thermodynamic equilibrium state of temperature T. If, however, the fluid is in the local-equilibrium steady state 
and can be characterized by the time-dependent temperature T = T{t), then ^ cannot be rigorously justified. We 
are optimistic and believe that it can be used as a first approximation to an exact (but non-existing) theory. In 
section 2, we present a method how to treat such a system when temperature can fluctuate with T{t) represented by 
a stochastic process and how to obtain the velocity distribution. In section 3, we consider a curtailed characteristic 
functional for which an evolution equation is determined by the infinitesimal generator of the stochastic process 
representing temperature fluctuations. In section 4, we consider a case of dichotomic temperature fluctuations and 
solve the evolution equation for the curtailed characteristic functional. In section 5, we discuss properties of the 
velocity probability density. In section 6, we analyze statistical moments of the velocity. 

II. TEMPERATURE FLUCTUATIONS: CHARACTERISTIC FUNCTIONAL 

Now, let us assume that temperature of the fluid fluctuates around its mean value Tq, 

T = T{t) = To + vit). (8) 

The zero-mean stationary stochastic process r]{t) describes temperature fluctuations and is independent of the stochas- 
tic process £,{t) describing thermal noise (interaction with surroundings). The formal restriction on this process follows 
from the condition T{t) > and its phase space Y is 

77(t)ey=(-To,oo). (9) 

The velocity probability distribution P{v,t) can be obtained from the relation Q, in which the initial transition 
probability density p{v, i|wo, 0) is expressed as 

p{v,t\vo,0) = :^ f cLue-''^^C,iu;,t;vo), (10) 

where the conditional characteristic function Cy{uj, i; vq) of the velocity is defined by 

a(c^,i;«o) = (e*""<*)). (11) 

The velocity v{t) is a solution of Eq. (Q with time-dependent temperature T = T{t), namely. 



v{t) = vq exp {—^t/m) + yj2k^i/m? I dsexp [— 7(t — s)/to] y/T{s) f(s). 







(12) 



where vq is the initial velocity of the Brownian particle. Inserting the above equation into (|ll|l yields 

Cy{uj,t;vo) = exp (iujvoe'^*^™] C{uj,t), 



where 



C{Lj,t) = ( exp 



Jo 



(13) 



(14) 



i,r, 



The subscripts ^ and ry denote average over all realizations of thermal noise £,{t) and temperature fluctuations ri{t), 
respectively. The averaging over the Gaussian white noise ^{t) can be performed leading to the expression 



C(w,i)=exp -{kTo/2m)uj^ (l-e-^''^^"') $^(w,i) 



where 



$^(Lj,t) = (exp 



-{'fk/m )uj e 



2V,,2«-27t/m 



ds e^'^/'^'Tjis) 



(15) 



(16) 



is the characteristic functional of the stochastic process T](t). In this approach, the velocity probability of the Brownian 
particle is determined by the characteristic functional of temperature fluctuations. The explicit form of this functional 
will be obtained by the method of the so-called 'curtailed' characteristic functional. 



III. CURTAILED CHARACTERISTIC FUNCTIONAL 

In order to calculate the functional IjlGjI we proceed in the following way ,1^J . For fixed time t — i we define [l^ 

n = (7fc/m')a;2e-2^*/'". (17) 

Let us introduce the auxiliary functional 



*[77;f^,i] = (exp 



~n I ds e'^^'I'^r^is) 
'o 



Then the relation 



holds. 



%{Lu,t)^^[f];n,i^t] 



The curtailed characteristic functional corresponding to (|18|) is defined as \lt 



V{y,t) = (S{Tj{t),y)exp 







(18) 



(19) 



(20) 



where y & Y takes values from the phase space of the stochastic process ri{t) and 6{ri{t),y) is the Kronecker delta 
when ri{t) is a discrete process or the Dirac delta for the continuous stochastic processes ri{t). The relation between 
these two functions is the following 



%ito,t)= / V{y,t)dy 



(21) 

where the integration for the continuous (or summation for discrete ) process is over the phase space Y . We intro- 
duce curtailed characteristic functional because for it, in contrary to l|16|) . an evolution equation is known. In the 
abbreviated notation, it has the form [13 



^Viy, t) = LViy, t) ~ fie^^*/™ yV{y, t), 



(22) 



where L is an infinitesimal generator (a forward operator) of the stochastic process ri{t). If ri{t) is determined by an 
Ito stochastic equation, the infinitesimal generator is a differential operator which occurs in the forward Kolmogorov 
equation (i.e. in the Fokkcr-Planck equation). Now, the problem reduces to solving the evolution equation (|22|l which, 
in dependence of r/(i), can be a single or a set of differential or integro-differential equations. Below, we present an 
example which can be solved exactly. 



IV. DICHOTOMIC FLUCTUATIONS 

Here, we consider a caricature of temperature fluctuations, i.e. a discrete, two-state model [l|. An extension to a 
many-state or continuous model of fluctuations is in principle possible jig . However, physics should be similar but 
mathematics would be much more complicated because of difficulties in solving the evolution equation (|22|l . So, we 
represent temperature fluctuations by dichotomic noise [iTJ 

Tj{t)^{~a,b}, 0<a<To, b>0. (23) 

Transition probabilities per unit time from one state to the other are given by the relations 

Pr{-a^b)=fi=l/Ta, Pr{b^-a) = X=l/n, (24) 

where Tq and T(, are mean waiting times in states —a and 6, respectively. We assume that 

bfi = aX. (25) 

Then the process is stationary and the probabilities 

Prirjit) ^ ~a) = ^ = ^ Pr^t) = b) = ^ = ^. (26) 

fi + A a + fi + A a + 

The first two moments read 

(77(t)}-0, {v{tMs))=abcxpi-\t-s\/T,), (27) 

where the correlation time Tc is given by the formula 1/tc — ^j, + X. 
The relation (|21(l takes the form 

<i>n{uj,t) = V{-a,t) + V{b,t) (28) 

and the explicit form of (|22|l reads 

|-V^(-a, t) = -^IV{-a, t) + XV{b, t) + ne^'^'/'^aV{-a, t), (29) 

^V{b, t) = ^lV{-a, t) - XV{b, t) - f]e2T*/"6F(6, t). (30) 

The initial conditions follow from H2U|I and read (cf. (|26|l ) 

F(-a,0) = (%(i),-a)) = ^, (31) 

V{b,Q)^{5{^{t)M^-^- (32) 

A* -I- A 

Let us define a new time variable 

r = T[t) = fle^-^*/™ (33) 

and two new functions V{—a,T) and y(6, t) via the relations 

V{-a,t)^V{-a,T{t)), (34) 

V{b,t)^V{b,T{t)). (35) 

Then Eqs (|30|l and (|29|l can be transformed to the form 

^^F(-a, r) = -^J.V{-a, r) + XV {b, r) + TaV{-a, r), (36) 

m OT 



5 

-^ — V(b, t) = iiVi-a, t) - \V{h, r) - rhVih, r) (37) 

m OT 



with the initial conditions dX T{t = 0) — Q., 



V{-a,n) = ^, (38) 



Vib,n)^^. (39) 

We define two new functions in the following way 

FiT)^V{-a,T) + V{b,T), (40) 

G{t) ^bV{b,T)-aVi~a,T). (41) 

Then from Eqs ^^ and ^7^ one gets 

^Fir) = -G{t), (42) 

TO 

— G(t) + U + \-T{a- b))G(T) + Ta5F(T) = 0, (43) 

TO 

where the dot denotes a derivative with respect to the argument. The initial conditions follow from (|40|l - (|43|l and 
take the form 

F(f7) = 1, G(fi) = 0. (44) 

The function F{t) is crucial because the characteristic functional (|16|) is related to it in a simple way. Indeed, 

^n{,uj,t)=F{T) for r = f^e^^*/" and ^ ^ {-ik/m^)uj'^e-'^^'''^ . (45) 

From the above system of two coupled differential equations (|42|l and (|43|l . we obtain a closed differential equation 
for the function F{t) only. It has the form 

4^, ^ m r , ^■, • . X m?abT , ^ , ^ 

tF{t) + — [^ + a + r(5 - a)] F{t) - -^:^F{t) = (46) 

with the initial conditions 

F{VL) = 1, F{VL) = 0. (47) 

It belongs to a class of hypergeometric equations. Its solution is the function [l3| 

F{t) = e-'"''-/^^ {Ci(17) $[a, /3, x(t)] + C^m *[«, /?, x(r)]} , (48) 

where $ and ^ stand for the confluent hypergeometric Kumnier and Tricomi functions, respectively 19]. The 
parameters 

mb a + X mb ,,„, 

a = — — — (49) 

2j a + b 27rc(a + 6) ' ^ ^ 

TO(^ + A)^^r^ (50) 

27 27re ^ ^ 



and 



, , m(a + b) , , 

XT = ^ r. 51 

27 



The constants Ciifl) and C2{^) are determined from the conditions (|47|l and read 

(a + ?))r(/3) 



and 



27 



TofO^ - ^^(") v(-0'|/3 p-maO/27 

X (ci>[a + l,/3+l,x(r!)]-$[a,/3,x(f^)]), 



(52) 



(53) 



where r(z) is the Euler Gamma function. In this way we found the function F{t) and via the expressions in 1451) we 
can find the characteristic functional $,,(1^, t). 



V. PROBABILITY DISTRIBUTION 

The probabihty distribution is obtained from Eqs 11U() - (|16|I and the relations H45|) . It has the form 
p(«,i|t;n,0) = -^ Hdoje^'^^ 



X exp 



The explicit form of the characteristic functional is 



(54) 



%{uj,t) 



br(a) 



^ — bkuj jlin 



[a + 5)r(/3) 
X i $ [a, /?, Ato"^] i * \a, (3, Auj'^e'^^*' 






27T-C 



+* [a, 13, Auj'^] $ [a + 1, /? + 1, Acj2e-27t/™l _ $ L^ ^^ Aw2e-27t/m 



(55) 



where the constant A ~ k{a + b)/2m. By use of l@J, the distribution p{v, t\vo, 0) allows to determine evolution of the 
one-dimensional velocity probability density P{v,t) for an arbitrary initial state defined by the distribution P(w,0) 
and analyze relaxation of the system to the stationary state. 



A. Stationary distribution 

The stationary velocity distribution function Pst {v) does not depend on the initial distribution P{v, 0) . It is obtained 
from Q and H54|l performing the long time limit, t ^ 00. We use the relations [l9| 



and 



*[a,/?,z] 



lim ^[a,l3, z] — I 



T{l-p) , r(/3-l) i_ 



r(a-/3+l) r(a) 



(56) 



(57) 



which represents the leading terms for small z = Auj'^e 27*/™ ^ \ when t ^ 00. Then the stationary distribution 
takes the form 



1 /""^ 



m a + b 



2^Tc{a + b) ' 27rc ' 2 m 



ku? 



(58) 



It is not a Maxwell distribution and we can conclude that the stationary state is not an equilibrium state: It is a 
nonequilibrium state. In the case of absent of temperature fluctuations, i.e. when a — and next 6 = 0, then 



$ 



2 rc7' 2jTc 
and the Maxwell distribution (O is obtained for T — Tq. 



,0 



= 1 



(59) 



B. Limiting cases of short and long correlation time 

The correlation time Tc of temperature fluctuations is deflned below Eq. H27|l . For very fast fluctuations when the 
correlation time is short. 



hm $ 



b m 



2jTa{a + b)'2jTc' 2 



m a + b , 9 



^bk cj /2m 



(60) 



and H58fl reduces to the Maxwell distribution (TJ with temperature T = Tq. The short correlation time limit can be 

achieved when (cf. (25)): 

(i) /i — *■ c» and 5 — > but bfi = aX = const. 

(ii) A ^ oo and a — )■ but aX = bfi = const. 

(iii) 6 — > (X) and A — *■ oo but b/X = const, (the latter corresponds to the Poisson white shot noise) 

(iv) ^ — > oo and A ^ oo but /i/A — const. 

In these limits, the system is not able to react to very fast fluctuations and effectively it feels the mean temperature 

T = To. 

The opposite limit is the adiabatic limit when fluctuations are slow and the correlation time is very long, Tc —>■ oo. 
The Kummer function takes the form 



Hm $ 



m a + b 



2jTc{a + b) ' 2jTc ' 2 m 
and Eq. (|58|l reduces to the function 

b 



kuj' 



Ja+b)kui^ /2m 



(61) 



P.t{v) 



a + by 2TTk{To — a) 
a 



a + b\l 27rk{To + b) 



exp 



exp 



2k{To - a) 



2k{To + b) 



(62) 



It is a linear combinations of two Maxwellian distributions for two temperatures Tq — a iTq + b and with the weights 
given by Eqs (|26l) . The long correlation time limit can be achieved when /i, A ^ and the mean residence times 
in the states —a and b tend to inflnity, Tq, r^ ~> oo. The adiabatic limit for other models of inverse temperature 
fluctuations has been considered in Q^ . 




FIG. 1: The stationary velocity distribution Pstiv) for three values of the correlation time Tc of temperature fluctuations: Tc = (dotted 
line, the Maxwell distribution), Tc = 5 (dashed line) and Tc ^ oo (solid line). 



VI. DISCUSSION 

Applying a method of the curtailed characteristic functional we obtained a time-dependent probability distribution 
of velocity of the Brownian particle moving in medium in which temperature fluctuates. We considered temperature 
fluctuations to be as simple as possible, i.e., a two-state stationary Markovian process r]{t). Nevertheless, the problem 
is formulated for an arbitrary Markovian stochastic process i]{t) because what we need is the infinitesimal generator 
L of the process 'ri{t), see Eq. (|22|) . We note that the correlation function of the force F{t) ~ y/2jkT{t) ^(i) in Eq. 
(Q has the form 

< F{t)F{s) >= 2-fkToS{t ~ s), (63) 

independently of statistics of fluctuations r]{t) and has the same form as in the case without temperature fluctuations. 
It resembles the dissipation-fluctuation relation. However, we showed that the stationary state is a nonequilibrium 
state. We can ask how far the system is from equilibrium. To this aim, let us analyze statistical moments of the 
velocity. From Eq. I|58() it follows that the stationary characteristic function Cu(a;) of the velocity is 



C,H=e-(^°+'')'="'/2'"$ 



m a + b , n 



2'^Tc{a + h) ' 27TC ' 2 TO 



' OJ 



(64) 



The statistical moments (w"), {n = 1, 2, 3, ...) can be obtained from the relation < (v") >= i^ cP Cy{uf) / duo'^\^=[). The 
odd order moments are equal to zero. The second moment 

(i;2) = fcTo/2. (65) 

It does not depend on the statistics of temperature fluctuations and is the same as for the Maxwell distribution! The 
forth order moment measure a deviation from the Maxwell distribution. We use it to calculate the kurtosis, 

where t^ = TO/7 is the relaxation time of the velocity in the deterministic case (cf. Eq. H12() when ^(i) = 0) and Tc 
is the correlation time of fluctuations (see Eq. l(23). For the equilibrium state, i.e. for the Maxwell distribution, the 
kurtosis is zero. In the case considered, the kurtosis is always positive and it means that Pst{v) is more peaked than 
the Maxwell distribution. It is an increasing function of the variance < 'rf'{t) >= ab and the correlation time Tc of 
temperature fluctuations. As a function of the correlation time, it grows from zero for Tc = to the maximal value 
ab/TQ when Tc — » 00. Generally, all even order moments are greater than for the Maxwell distribution. 

One can determine the moments for the position of the Brownian particle. E.g., for long times, t ^ r^, the mean 
squared displacement 

{x^{t)) r^2Dt, (67) 

where the diffusion coefficient D — kTo/"f is the same as in the case without temperature fluctuations. It means that 
for long time, the process in the position space is the standard normal diffusion with the same diffusion constant. 
We can conclude that the flrst two moments of position and of velocity are the same in both cases: without and with 
temperature fluctuations. So, when we measure only flrst two moments, we can not distinguish these two states. We 
emphasize that it does not depend on the model of temperature fluctuations. 
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